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ABSTRACT 



In order to test the unification scheme for double- 
lobed radio sources, the far-infrared properties of 
matched samples of radio-galaxies and radio-loud 
quasars were studied using ISOPHOT. The quasar 
data were complemented with nearly simultaneous 
submillimeter and centimeter radio data. The re- 
sults show that quasars are generally brighter than 
radio-galaxies in their far-infrared output and that 
beamed non-thermal radiation must be excluded as 
source for this infrared excess. However, non-thermal 
flares or optically thick tori can still reconcile these 
results with the unification scheme. 

Key words: active galactic nuclei; unification; spec- 
tral energy distributions 



1. INTRODUCTION 



The unified t heorv for radio-loud ac tive galaxies 
dBarthel, 1989| ; |Urry fc Padovani, 1995| ) predicts the 
existence ot a dusty torus surrounding the central 
engine. Models for the reprocessing of the hard non- 
thermal nuclear radiation by this torus show that at 
wavelengths longwa,rd of BO/zm the torus will be op- 
ticall y thin ( Pier fc Krolik, 1992 ; Granato fc Danese, 
1994 ) . However, in the IRAS observations quasars 
are two t o four times brighter at 60 fiui than radio 
galaxies ( Heckman et al., 1994 ). In some special ge- 
ometries of the torus, models can account for this 
infrared excess. Also the torus could still be op- 
tically thick at 60 /im, or there could be a signifi- 
cant amount of angle dependent non-thermal con- 
tamination in quasars. Non-thermal beamed radia- 
tion has been proposed as a possible solution (Heset 
al ., 1995| ). 



In order to discover the true cause of the infrared 
excess in quasars, observations are needed at longer 
wavelengths. All current torus models predict that 
beyond 80 /im the thermal emission will be isotropic. 
According to unified schemes, the output of quasars 
should then equal that of matched radio-galaxies at 
these long far-infrared wavelengths. ISOPHOT ob- 
servations of quasars and radio-galaxies have been 
obtained to test this hypothesis. To quantify any 
possible non-thermal contamination in quasars, the 



latter were in addition observed with SCUBA at the 
JCMT, and the NRAO Very Large Array (VLA). 



2. SAMPLE SELECTION AND OBSERVATIONS 



The target sample is selected using several crite- 
ria. First of all the objects must have the edge- 
brightened, double-lobed radio str ucture typical for 
Fanaroff & Riley class II sources (Fanaroff & Riley, 



1974). Six matched pairs were constructed, contain- 



mg one quasar and radio-galaxy each, having com- 
parable redshift, radio-lobe power and angular size. 
Lobe radio power is a reasonably good measure of the 
active nucleus power, and is furthermore independent 
of the orientation of the objects. 

Observations were carried out with ISOPHOT at 
three wavelengths, during the period October 1996 
to April 1998. We made raster- mode mini-maps with 
ClOO and C200 (AOTs P3 and P22). The rasters 
have a size of 3 x 3 (ClOO) and 4 x 2 (C200) in 
YxZ direction. The filters used were the 60 /im, 
90 /um (ClOO) and 160 /im (C200). Observations of 
only four pairs were completed during the lifetime of 
ISO; for the other two pairs only the radio-galaxy was 
observed. For the quasars these data were comple- 
mented with nearly simultaneous VLA A-array data 
at 6, 2, 1.2 and 0.7cm (April 17, 1998) and JCMT 
SCUBA data at 2mm and 850 ^im (May 28, 1997), 
addressing their non-thermal radio core strengths. 
These observations were made in order to detect pos- 
sible non-thermal contamination of the far-infrared 
emission. Since the effect of non-thermal emissio n 
in radio-galaxies is negligible ( Hoekstraet al., 1997 ), 
these objects were not observed with the VLA and 
SCUBA. 



3. DATA REDUCTION 



3.1. ISO Data Reduction 



All ISO data have been processed with OLP version 
6.11 and raw data have been reduced using the Phot 
Interactive Analysis tool (PIA) version 7.2(e). The 
steps from ERD to AAP level were all done using 
the default procedures. Deglitching at ERD level was 



Table 1. Flux densities in mjy as plotted in Figure 1 (table added after submission to proceedings) . For 3C 351 
the flux density at 6 cm is from our VLA observation, while in Figure 1 we plotted the datapoint from Bridle et 
al. (1994). 
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done with the two-treshold method, setting sigmas of 
3.0 (flagging) and 0.5 (reacceptance). At SRD level 
reset interval correction was applied, as well as dark 
current subtraction per raster point. Afterwards the 
data were deglitched again with a sigma of 2.5 (source 
measurements) or 0.5 (FCS measurements). At SCP 
level no significant corrections were applied. The sig- 
nal of the FCS was in most cases outside the hard 
extrapolation limits of PIA, therefore it was decided 
to use only default responsivities for all sources at all 
wavelengths. 



3.2. ISO Flux Determination 



The final flux was determined using an IDL program 
called ana_niini_map, developed by the PHOT sup- 
port center in Heidelberg. This program uses each 
pixel of the detector as a scan across the source. 
Thus no map is made, instead it determines the back- 
ground flux from the two off-source raster positions 
immediately before and after the on-source raster po- 
sition. This off-source flux is compared to the on- 
source flux. If a signal is present, all pixels will show 
a positive difference between on- and off-source po- 
sitions. However, due to the faint flux levels, noise 
caused a lot of false results for some pixels. To avoid 
using noisy data, an IDL tool was developed to deter- 
mine the noise on the pixels used by ana_mini_map. 
Pixels with extremely high noise on either one of the 
off-source positions or the on-source position were ex- 
cluded from the flux determination. This produced a 
reliable flux for all the detections (errors are generally 
less than 10 per cent). 

Tests were conducted by comparing the results to 
DIRBE background measurements and IRAS detec- 
tions of our sources. Both DIRBE and IRAS repro- 
duced our results within a 10 per cent error limit. 
However, due to the large calibration errors for faint 
sources, the fluxes presented h ere are to be take n 
with errors of about 30 per cent (Lemke et al., 1996). 
In order to present the results as clearly as possible, 
only the statistical errors are given by the error bars 
in the plots, thus providing better insight into the 
real meaning of the detection. 



3.3. VLA And SCUBA Data Reduction 

The VLA data were reduced using the AIPS version 
of 15 April 1998. Calibration was done using the 
standard flux calibrators and nearby phase calibra- 
tors. After calibration the maps were cleaned us- 
ing CLEAN and the core flux density values were 
determined using the MAXFIT procedure. 3C351 
caused some problems at 6 cm and 2 cm, because of 
the strong hotspots in this source. The exact flux of 
the core in these band s has not vet been determined. 
The 6cm core flux of Bridle et al., 1994 was plotted 
instead, since no strong variations have been reported 
of this source at cm wavelengths. For all sources the 
errors are smaller than the symbols in the plots. 

The SCUBA data were reduced using the SCUBA 
Data Reduction Facility (SURF) and Kernel Appli- 
cation Package (KAPPA). Not all sources were ob- 
served at all bands. Only sources with an observation 
have a SCUBA upperlimit, and only 3C 334 has been 
marginally detected (~ 2cr), be it in both bands. 



4. RESULTS 



The resulting broad band spectra are plotted for 
each pair in Figure 1. It is immediately clear that 
in most cases the quasar is still brighter than the 
radio-galaxy, even up to the longest ISOPHOT wave- 
lengths, with only one exception. The quasars are 
detected at all ISOPHOT bands, whereas in only 
one out of four observations the radio-galaxy was de- 
tected. From the additional VLA and SCUBA data it 
is evident that the quasar cores do not contribute sig- 
nificantly to the far-infrared emission observed with 
ISO, since the sub- mm points and VLA 7mm data 
would have to be much higher in that case. 

From the simple black-body fits to the ISO data it 
is also clear that two of the four radio-galaxies need 
a lower temperature than the quasars. This could 
imply that in radio-galaxies we observe a cooler dust 
component, which is in agreement with the unifica- 
tion theory, assuming a dust distribution with a tem- 
perature gradient. In one case the black-body tem- 
perature could not be determined, because there is 
no detection of the radio-galaxy at any wavelength 
(3C19). In the final case the far- infrared spectra of 
the quasar and radio-galaxy are identical, which is 
surprisingly enough the pair with compact sources. 
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Figure 1. The complete spectra for the four pairs observed with ISO. For the quasars: plusses are VLA data, 
triangles are SCUBA data, red diamonds are ISO data. The solid line is an empirical fit to the radio-core spectra, 
the dash-dotted line is a black body fit to the ISO data. For the radio-galaxies: asterisks are ISO data and the 
dotted line is a black-body fit. 



5. DISCUSSION 



5.1. Non-thermal Contamination 



Radio-quiet and radio-loud quasars show remarkable 
similarities in their broad band spectra. The simi- 
lar infrared to submillimeter slope and far-infrared 
turnover frequency point to an identical source of 
infrared emission. Since radio-quiet quasars have 
no significant non-thermal component, the inf rared 
emission must be thermal ( Sanders et al., 1989 ) . As- 
suming the models for thermal torus emission show 
the correct trends, an aspect angle dependent non- 
thermal contamination must be present in radio-loud 
quasars to reconcile the observed infrared excess in 
these objects with the unification theory. There 
are two possibilities: continuous beaming, or flaring, 
which requires a rapidly varying beamed component. 



Continuous beami ng is in most cases no t more than a 
30 per cent effect (Hockstra ct al., 1997) and has been 
solidly ru led out for causing the exc ess in at least two 
quasars (van Bemmel et al., 1998). Here, only one 



quasar is detected at 7 mm and the detections with 
SCUBA are marginal at best. In our results continu- 
ous non-thermal contamination of the far-infrared is 
negligible, as can also be seen from the fits to the ra- 
dio data in Figure 1. It is important to use only core 
fl-Ux densities, as this is the only source of beamed 
non-thermal radiation in active galaxies. When us- 
ing the total radio flux densities as done by Haas et 
al. (1998), the centimeter flux densities will be dom- 
inated by emission from the radio lobes. These have 
generally steep spectra and are not beamed, thus 
they do not contribute significantly to the infrared 
emission in these sources. 

A varying component can easily be missed. Strongly 
beamed radiation fr om non-thermal fla res has been 
observed in blazars ( Brown et al., 1989| ); such flares 
are both very strong and short-lived, with typical 
timescales of several days. They could be observed 
in quasars due to the small aspect angle of these ob- 
jects. In previous work two out of three quasars 
showed high 3 mm fluxes, which could not be fit- 
ted wit h a radio core model and were suggestive of 
flaring (van Bemmel et al., 1998). In radio-galaxies 
beamed flare radiation is much less likely to be ob- 



served, because of the larger aspect angles in these 
sources. Strong flares can cause large differences in 
far- infrared output, whilst being missed in the non- 
thermal core spectra. 



5.2. Optically Thick Torus 



The persistence of the infrared-excess at longer wave- 
lengths can also be explained if the torus is optically 
thick up to 140 /im at least. Arguing in favour of 
this concept is the fact that a lower temperature is 
needed for the black-body fit in at least two of the 
radio-galaxies. Here maybe only the cooler edge of 
the torus is seen, while in quasars the hot inner parts 
dominate the infra-red output. This is in contradic- 
tion with the present models for dusty tori, but in 
broad agreement with predictions from the unified 
models. Since changing the models for the torus also 
has implications for the X-ray spectra observed in 
these sources, new modelling has to be done to study 
this possibility in more detail. 



5.3. Evolution Between Types 
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Evolution has b een proposed as a solution to the 
infrared-excess ( van Bemmel et al., 199q ). This 
would then imply that quasars are younger and still 
contain more dust in their host galaxy, whereas radio- 
galaxies are older and have processed their dust and 
gas into stars. A link with the ultra-luminous galax- 
ies has been made in this respect as these be ing very 
voung quasa rs who are in the phase of birth ( panders 
et al., i988| ). However, during a recent conference 
on this topic, most evidence showed that these so- 
called ULIRGs are in many cases driven by strong 
starbursts and not by a central active nucleus (Mon- 
sters or Babies: What powers ultra-luminous infrared 
galaxies?, Ringberg Castle, 1998). The fact that we 
find one pair in which the far-infrared spectra are 
identical, also argues against evolution between the 
types. Any evolutionary model is obviously in con- 
fiict with orientation based unification concepts. 



6. CONCLUSIONS 



Observations of matched pairs of quasars and radio- 
galaxies show that the previously reported infrared- 
excess is still present in quasars, up to 140 /xm rest- 
frame wavelength. This cannot be accounted for by 
current torus models, without causing a contradic- 
tion with unified schemes. Continuous beaming has 
to be excluded in our objects, but non-thermal flares 
can explain the excess, without necessary adjust- 
ments to either theory. Also the models for the torus 
have to be examined in more detail to explain ob- 
served temperature differences between quasars and 
radio-galaxies. This is currently under investigation 
and will be described in future publications. 

Future work will also include the search for flares in 
3C 334 with SCUBA and the investigation of orienta- 
tion dependence of optical emission lines in a sample 
of double-lobed active galactic nuclei. 
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